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Abstract 



We considered propagation of Ultra High Energy Cosmic Rays (UHECR) 
through the galaxy. We investigated models with sources of UHECR distributed 
in the same way as Cold Dark Matter (CDM) in a self-consistent way, taking 
into account both extra-galactic and Galactic contributions. Using a very sim- 
ple toy model of galactic magnetic field we showed that in the case of galactic 
origin of UHECRs the anisotropy can reach considerable values. In the case of 
extragalactic UHECRs origin the anisotropy appears to be non-vanishing only 
for electron and photon components due to synchrotron losses, but it hardly 
can be reassured. The reason is an extremely low flux of UHE electrons and a 
too low level of 7-ray anisotropy. 
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1 Introduction 



The observation of Ultra High Energy Cosmic Ray (UHECR) events with energies 
above 10 20 eV [II], Q is one of the most intriguing puzzles of modern physics. The 
point is that if UHECR come from large distances their energy dramatically decreases 
due to interactions with background fields. For example, nucleons of energies above 
~ 4 x 10 19 eV start to interact with the Cosmic Microwave Background (CMB) 
and lose their energy due to pion photoproduction (the so-called Greisen-Zatsepin- 
Kuzmin (GZK) effect |J). As a result, nucleons can not travel to an observer from 
distances larger than Rgzk — 50 Mpc. Heavy nuclei travel even much less distances 
of about few Mpc 0] because they photodisintegrated by CMB. UHE electrons lose 
their energy due to synchrotron radiation in extragalactic magnetic fields. From other 
side, it is hard to accelerate protons, electrons and heavy nuclei up to such energies 
even in the most powerful astrophysical objects such as radio galaxies and active 
galactic nuclei, but even if it is possible, usually such objects are situated far away 
from our galaxy on distances larger than Rgzk- This means that energy of a UHE 
particle will drop below the GZK threshold unless its energy is larger than observable 
one on several orders of magnitude. 

It was suggested several ways to avoid this difficulty, but all of them require 
either the introduction of a new physics well beyond the Standard Model or extreme 
astrophysics. 

One possible way to avoid GZK cutoff is to introduce a new particle, which inter- 
actions with CMB are suppressed in compare with the proton of the same energy. The 
possible examples of such particles are light long-leaving gluino-containing baryons 
U or axion-like particles 0. 

Another possibility is to stand with Standard Model particles but change physics 
laws. For example, protons and neutrons still can be UHECR with energies larger 
then cutoff if one slightly violate Lorentz invariance ||. 

One more way avoid a GZK cutoff is to produce UHE protons and photons inside 
of GZK volume from hidden sources. Photons with extremely high energies E > 10 22 
eV have attenuation length larger then 100 Mpc, if Extragalactic Magnetic Fields 
(EGMFs) are small -Begmf < 10~ n G. Such photons can be the sources of the sec- 
ondary UHE photons inside of GZK sphere |J. Another possible source of both UHE 
protons and photons are neutrinos with E > 10 22 eV, which could interact with back- 
ground neutrinos and produce secondary particles through Z° resonance [ID]. Both 
those scenarios require extreme astrophysical sources which should accelerate primary 
photons up to energies E > 10 23 eV. 

The other way to explain UHE events above the GZK cutoff is to create UHE 
protons, leptons and photons in decays of supermassive X-particles inside of the GZK 
sphere. This is one of the so-called 'top-down' (TD) scenarios. Depending upon the 
particular scenario, decaying X-particles could populate the Galaxy or can decay on 
cosmological distances In the latter case the important role is played by their 
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decay secondaries' energy losses during propagation through the space. The sources 
of the massive X-particles themselves might be topological defects such as cosmic 
strings or magnetic monopoles that could be produced in the early Universe during 
symmetry-breaking phase transitions envisaged in Grand Unified Theories |T3]. In 



an inflationary early Universe, the relevant topological defects could be formed at a 
phase transition at the end of inflation. Other possibility is to produce superheavy 
long-living particles thermally during reheating epoch of the Universe (13| or from 
vacuum fluctuations during or after inflation |fTI| . For a review of TD models see, for 
example fTSfl . 

Important information about cosmic ray sources give the measurement of the 
UHECR arrival directions. All experiments see uniform distribution of the UHECR 
over the sky. This leads to the conclusion that most probably UHECR come from 



the extragalactic sources. Besides, AGASA experiment see small scale clusters |JT6 



Those clusters can come from point-like sources |I7|. Moreover, the correlations of 



UHECR with BL Lacertae was found . If those small scale clusters and especially 
correlations with BL Lacs will confirmed in the future, TD models will be excluded. 

The main goal of this work is a self-consistent consideration of the galactic and 
extragalactic contributions to spectra of UHECR from TD models. There are two 
major possibilities: if sources of UHECR are present in our Galaxy or not. In the 
first case simple estimates show that contribution of the Galaxy component is about 
100 times larger [|19|] than extragalactic component and hence one might take into 
account the Galactic contribution only. Our numerical calculations, based on detailed 
modeling of UHECRs propagation process give even more strict limitations on the 
ratio of the extragalactic component of UHECRs to the galactic one. 

The main observational signature allowing to distinguish this kind of models is a 
quite significant anisotropy in all the UHE particles spectra which is in fact simply 
due to the non-central position of the Sun in the Galaxy at the distance 8.5 kpc 
from the Galaxy center. In what follows we present our results for various models of 
UHECR source distributions in the Galaxy. 

In the case there is no sources of UHECRs in the Galaxy, the anisotropy in all 
particle spectra is vanishing, except for photons and electrons, for which a small 
anisotropy in the directions to the Galactic center and opposite one still exists due 
to synchrotron radiation. 



2 Propagation through the Galaxy 

A numerical codes has been developed recently for calculation of spectra of UHECRs 
in the case of isotropic and homogeneous source and magnetic field distributions 



D| |20, [Hp . Our code provides a rather reliable information on possible extragalactic 
sources of UHECRs. Unfortunately, as was pointed earlier |T| , the method developed 
in above papers may not be used directly if one is interested in sources within our 
Galaxy, moreover the galactic magnetic fields may not be neglected even for the 
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extragalactic component of UHECRs in the case of light charged particles (electrons, 
positrons). Electrons loose their energy rapidly enough due to synchrotron radiation 
in the galactic magnetic field: 

dE 4 a 2 B 2 E 2 



dt 9 mt 



e 

-6/ 




where B ~ 10 G is the galactic magnetic field, and E is the UHE electron energy. 
Then the energy loss length is 

= 3 - 8x (Toivf (lo^f «* c - (2) 

Deflection is another important factor when dealing with the propagation problem 
in a non-anisotropic case. The straight line propagation (SLP) approximation which 
treats the motion of CR particles in one dimension fails if the effect of the deflection 
becomes large. The gyroradius of a charged particle with charge qe and momentum 
p (energy E) is given by 

P E 1 / E \ ( B ± \ _1 

R ° = i^T " ^b! " 11 x ^lib^vj Iio^gJ kpc ' (3) 

where B± is the field component perpendicular to the particle's momentum. 

However, if one is interested only in the spectra of particles with energies larger 
than 10 19 eV, i.e. particles which weakly (up to 10 degrees) deflect in the galactic 
magnetic fields B ~ 10 _6 G on the distances lOkPc, one can use a straight-line prop- 
agation approximation. This means that one can use the same kinematic equations 
as in JH], ^] but written for the differential fluxes instead of concentrations. 

j t 3*(K, t) = ~3a(E a , t)a a (E a , t) + (4) 



£ J (3 ca (E c , E a , t)j c (E c , t)dE c + 



Q a (Ee,t) 

47T 



where j a (E a ,t) is a differential flux of particles a. 

Integration in Eq. (f|) is performed on the chosen straight trajectory within the 
Galaxy. The initial values of j a {E) may be taken from the previous task (extragalactic 
component). A some another difference is in a time dependence of external source 
and magnetic field terms. Their values are determined by current position within 
the trajectory. Further, for simplicity we will assume spherical symmetry, that is the 
source concentration and the magnetic field strength are taken to be functions only 
of distance from the center of the Galaxy r. 

B(t) = B'(r(t)), Q(t) = Q'(r(t)), (5) 
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Figure 1: Propagation in the Galaxy. Points C and E show locations of the galactic 
center and Earth respectively. 

As a toy model of galactic magnetic fields we will use the one decreasing exponentially 
with r: 

B\r) = B exp(—r/r B ). (6) 

In the Fig{l] we plot point E, corresponding to the Earth location at a distance a 
from the galactic center C. Here R is an effective Galaxy radius, the distance at which 
the source concentration and magnetic field strength become equal to extragalactic 
values. Let us draw an arbitrary trajectory crossing point E (horizontal line in the 
Fig.|l|), a being the angle between the direction to the center and the trajectory and 
t being the length (time) parameter on the trajectory (point t = is the nearest one 
to the center of Galaxy). 



r(t) = V a 2 sin 2 a + t 2 , (7) 
For the flux coming from "center" limits of integration are given by 



tmin = -\AR 2 — a 2 sin 2 a, t max — Va 2 — a 2 sin 2 a; (8) 
for the flux from the opposite side 

t min = V a 2 — a 2 sin 2 a, t max = V R 2 — a 2 sin 2 a. (9) 
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3 Source distributions 



As we have already pointed out in Introduction, the absence or presence of sources of 
UHECR in the Galaxy lead to quite different UHE particle spectra on the Earth. In 
the first case only electron and photon spectra are slightly anisotropic, while in the 
latter case an anisotropy is there for all particle spectra. 

We considered the same models of source distributions as in 111911 . First is the 



isothermal halo model 22 



™( r ) 00 / 2 ! 2 N > ( 10 ) 

[r A c + r z ) 

where r c is some core size. A model with the distribution ( |T0"D may describe the 
collisionful CDM case. The second model has the distribution of Ref.[Q, which 
corresponds to the opposite limiting case of collisionless CDM 

n(r) oc = . (11) 

(r 2 + r 2 )(R + r) 2 

We will take always r c = 2kpc and R = lOOkpc. 

Since the distribution of CDM in the Galaxy might contain most probably both 
collisionful and collisionless components, the results for anisotropy in a realistic case 
is probably somewhere in between of two extreme cases we considering here. 

For the ratio of mean galactic and overall X-particle concentrations we used value 

< UX - 1 > = 1 = 10 5 . (12) 
< n x > n gal V ga i 

4 Results 

The results of our numerical calculations of propagation of UHECR through the 
Galaxy for the model of long-living X-particles |T3] are presented on Figs. 2-6. As 



an example, we used the following parameters of the model: we took mass of a 
superheavy particle m>x = 10 23 eV = 10 1A GeV. We took as the main decay channel 
the X — > qq one. We assumed that X-particles give the main contribution in the Cold 
Dark Matter (CDM) with Vt x = 0.3 (assuming fi A = 0.7 and H = 70 km s^Mpc -1 ). 
We normalized the resulting UHECRs flux to the observed one (AGASA data), and 
calculated required tx using the fact that fluxes ji ~ t x . Extragalactic magnetic 
fields were taken to be equal to B out = 10~ 12 G. We used usual QCD-motivated 



fragmentation functions (see Eq. (52) in [15|) 



We defined anisotropy in spectra of UHECR in the following simple way: 

< „ ^center ^border /i o\ 

= 1 Tl ' ^ ' 

•-'center > '-'border 
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Figure 2: Spectra in the long-living X-particle model with ttx — 0.3, = 0.7, 
H = 70 km s~ 1 Mpc~ 1 , t x = 8.6 x 10 7 t[/ (% is the age of Universe) under constant 
EGMF of strength B = 10 _12 G, assuming decay mode X —>■ qq. This spectra is taken 
near the boundary of the Galaxy. The observed spectrum at energies more than 10 19 
eV will be similar, but without electrons, which lose all the energy due to synchrotron 
radiation. For this reason we don't show here the electron spectrum (in fact it nearly 
coincides with photon one at high energies). 

where J C enter and Jborder are the differential fluxes of UHECR from the direction on 
the Galactic center and from opposite direction. 

We considered two essentially different cases of the X-particle distribution. In the 
first case we took uniform distribution of the X-particles in the Universe, but without 
X-particles in our Galaxy. In the second case we considered X-particle distribution 
similar to cold dark matter distribution in the galaxies. 

On the Fig.|2] we plot spectra of UHE protons and photons near the boundary 
of our Galaxy in case of extragalactic sources of UHECRs. Propagation through 
the Galaxy mostly leaves the picture at highest energies without changes. It is not 
so only for electrons, since they are strongly suppressed by galactic magnetic field. 
Lifetime of X-particle in this case is tx = 8.6 x 10 7 tjj (tjj is the age of Universe). 
From Fig.§ one can conclude that all highest energy cosmic rays are photons. This 
is the general result for all TD models. Moreover, for models with mj > few x 
10 22 protons hardly can contribute significantly in the region below cutoff, since an 
attempt to increase extragalactic magnetic field to fit simultaneously highest energy 
part of AGASA spectrum by photons and lower energy part by protons will lead to 
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Figure 3: Anisotropy in spectra of photons due to synchrotron radiation on galactic 
magnetic fields for the model of Fig.[| Galactic magnetic field is taken in the form 
(|6[), where B = 1CU 6 G, re = Akpc. There is no anisotropy in proton spectrum in 
this case, because protons does not lose the energy. 



overproducing of photons at EGRET region, due to synchrotron radiation. And so 
in this case one need to suppose that another kind of UHECR sources exists in order 
to explain UHECR below cutoff. 

The scale of spectra distortions produced during propagating through the Galaxy 
can be illustrated in terms of the flux asymmetry in directions to and from the Galaxy 
center Eq.(|T^). It is presented on the Fig.[5| The spectra in any directions are similar 
for all particles, except for photons: they are a little bit different due to influence 
of galactic magnetic fields. Namely, some part of photons, which was converted to 
electrons, was lost due to electrons synchrotron radiation in Galactic magnetic field. 
Because the path through the Galaxy center is longer then from opposite direction, 
but both of them much smaller then photon mean free path, there is small asymmetry 
in photon flux, as we see from Fig|J. 

At Fig. | we presented typical UHECRs spectra in the case of the galactic origin. 
As was expected, the UHECR flux in this case is determined by the source distribu- 
tion rather than by propagation processes. This is not the case only for electrons, 
because of the extremely small synchrotron energy loss length in magnetic fields of 
the Galaxy. The A-particle lifetime is tx = 1.5 x lO 11 ^. We checked that con- 
tribution of UHECR from all other galaxies is negligible in this case in compare to 
the contribution of the Milky Way. The approximate relation for the Galactic versus 
extragalactic fluxes is given by relation of corresponding X-particle lifetimes. So, we 
get F G aiaxy /-^extragalactic ~ 1.7 x 10 3 . As we see from Fig. |, the spectrum of UHECR 
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Figure 4: Spectra of UHECR coming from Galaxy center in the case of (^Tj) source 
distribution (similar to collisionless CDM distribution) in Galaxy, tx = 5 x lO 11 /^. 
The electron flux is not shown here, it is at least 100 times less then other fluxes 
at energies more then 10 19 eV. Flux of UHECR is dominated by photons, however 
another source of UHECR with energies below GZK cutoff is required. 
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Figure 5: Anisotropy in spectra of UHECR in directions to and from Galaxy center 
for the thermal source distribution in Galaxy (collisionfull CDM). 



dominated by photons at all energies. Protons contribute at most up 10% to the flux. 
This is simply due to fragmentation of quarks in X-particle decay. 

In the case of Galactic origin of UHECR the important parameter is anisotropy in 
particle spectra, defined in Eq. QT5|). We investigated two essential cases of collisionless 
Eq.flllD and collisionful Eq.([H]) dark matter. The asymmetry in particle spectra for 
those models is plotted in Fig.|5| and Fig.|| 

We conclude that anisotropy in spectra of nucleons, photons and neutrino is an 
important test for top-down models, in which sources are located in galaxies. In the 
realistic models of source distribution the anisotropy can rise up to 80-90%. The 
dominant contribution in this case is provided by a flux coming from our Galaxy. 
The extragalactic input in this case appears to be at least 1000 times smaller. The 
shape of the spectra at energies higher than few xlO 19 is determined mostly by 
source distribution. Deflection in galactic magnetic field is negligible and so the exact 
configuration of the fields is not important. The flux of UHECR will consist mostly 
from photons at highest energies. 

In the case of the extragalactic UHECRs origin there will remain some anisotropy 
in photon and electron spectra, however extremely low flux of electrons and low 
value of photon anisotropy (less than 1%) make the anisotropy hardly detectable. 
Nevertheless, it is probably worth to notice that in this case more exact models of 
galactic magnetic fields may change the final result somewhat. 
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Figure 6: Anisotropy in spectra of UHECR in directions to and from Galaxy center 
for the ( |TT| ) source distribution (similar to collisionless CDM distribution) in Galaxy. 
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